Goddard Space Flight Center 


Generalized Linear Covariance Analysis 

J. Russell Carpenter 
F. Landis Markley 


Motivation 




x(t) 

P G _, R,Q 

\ 


Order 

Reduction 

Likely 


s(t) = S (t)x(t) 



Error Budgets 


A A A 

P„-, R,Q 


Implementable 

Model 




d 


Generalized Linear Covariance Analysis 


2 


Carpenter and Markley 




Goddard Space Flight Center 


Some Background 


• Jazwinski (1970) - “solve-for/consider” state framework; 
assumes sequential filter 

• Gelb (1974) - “true/filter” state framework (“true state” is 
“difference” between filter state and truth state); assumes 
sequential filter 

• Maybeck (1979) - variation on Gelb; uses true (linear) state, 
rather than deviation from filter state 

• Markley, et al. (1988) - uses solve/for consider framework, 
explicitly models contributions from a priori, measurement 
noise, and process noise for batch and sequential 

• Tapley, et al. (2004) - uses solve-for/consider framework for 
batch and sequential 
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Present Work 
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• Update to Markley et al. (‘88, ’89) 

• Explicitly addresses sensitivity of solve-fors with respect 
to considers in context of ‘88 paper’s covariance partitions 

• Explicitly addresses “postdiction” in the batch framework 

• Applies method to integrated orbit/attitude problem 

• Includes some new ways to examine output 
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Context ... 
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Actual Trajectory: X(t) 
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Linear Deviations: 
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Context, Continued 
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Solve-for and Consider State Mapping 


Solve-For States “Consider” States 

s(t) = S (t)x(t), c(t) = C (t)x(t) 


Inverse Mapping: 

s (t) 


M(t) = 


C (<) 


, M-'(t) = [S(t), C(i)] 


x(t) = S(t)s(t) + C (t)c(t) 
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Covariance Partitions 
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Effect of a Priori Errors: AP a = SP 0 S T — P Q 


Effect of Measurement Noise: AP V = $P„S T — P v 
Effect of Process Noise: AP U , = SP U ,S T — P u , 

N.B.: A’s are not necessarily positive or negative 
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Auxiliary Notation 
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Solve-For Mapping of Dynamics and Measurements Partials 


^ss(ti,to) = S(ti)4>(ti,t 0 )S(t 0 ), H s (ti) = 


Process Noise Covariance Matrices 
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Sequential Filter Propagation 


Only the “a Partition” receives the a Priori Errors ... 
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... and only the “if Partition” receives the Process Noise 
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Sequential Filter Update 



The Gain is computed from the Formal Solve-for Covariance: 
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Only the “v Partition” receives the Measurement Noise 


Generalized Linear Covariance Analysis 


11 


Carpenter and Markley 


Batch Estimator without Process Noise 



We write the batch update in a form resembling that of the 
sequential filter to isolate the contributions of the a priori error and 
the measurement noise: 


Ki = [p*i + J (ij ) ^7" H s )4> ss ,*,)]“ 1 <I>I S («*,«.) HJ (**) ftr 1 



IT 


We associate the a priori error with an anchor time that is not 
required to precede the measurement batch. 
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Effect of Process Noise on Batch Estimator 
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The batch estimator does not model process noise, but 
if it is present, it will induce an error: 
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Generalized Process Noise Covariance 
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Process noise enters the innovations, and thus must be mapped 
to the anchor time: 


Qd (t* 5 tj ) 
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If the anchor time is between a pair of measurements, their 
process noise contributions do not overlap, so the expected 
contribution to the covariance is zero 
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Sensitivity to Individual a Priori Parameters 


• In this work, a “sensitivity” is a matrix of partial derivatives of 
some parameters of interest with respect to others (cf. Gelb, 
where “sensitivity” refers to a covariance matrix) 

• The “AP’s” give an indication of sensitivities of solve-for 
covariances to groups of errors (a priori , measurement noise, and 
process noise) 

• Often we’d also like to know the sensitivity of each solve-for at 
any point in time to each individual a priori solve-for or consider 
error 

Sensitivity of solve-fors with respect to a prior? s (sequential form): 

Z a (ti) = [I - Z a (t 0 ) = [S(i 0 ),C(i 0 )] 

Sensitivity of solve-fors with respect to a prior? s (batch form): 

Z a (t) = S(t)<l>(M.)[l - £LiS(^K;H(^(^M[S(* *J,C(M] 
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Example: Entry Problem 
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Altitude 


We want to 
understand the effect 
IMU parameters have 
on navigation 
performance 
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IMU Model 
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Accelerometer state : acceleration bias, scale 
factor, non-orthogonality 


Sac = [ b D(ac) F(a c ) ] x a 

= H a x a . 
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Variance Sandpile: F-Component of Inertial 
Position Error 
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Variance Sandpile: F-Component of Case- 
Fixed Gyro Angular Error 


Variance Sandpile 
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Solve-Fors at Final Time 

Gyro A SF Abias Vel Pos 


Sensitivity Mosaic 
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Logarithmic Sensitivity Mosaic (dB) 

A Priori State Index 
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Summary 
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Present work updates 
to Markley et al. 

Augments approach 
to sensitivity analysis 

Addresses 
“postdiction” in the 
batch framework ~ 

Applies method to 
integrated 
orbit/attitude problem 

Includes some new 
ways to examine 
output 
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IMU sensed acceleration error integrates into postion and velocity 
errors according to the INS model: 
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Example Problem Parameters 
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Simulation Parameter 

Value 

Units 

Gravitational Constant 

4.305 x 10 4 

km 3 /see 2 

Measurement Time Interval 

2 

sec 

Estimation Parameter 

Standard Deviation 

Units 

True Position Measurement Noise 

0.305 

m 

Formal Position Measurement Noise 

0.914 

m 

Initial Position Error 

30.5 

meters 

Initial Velocity Error 

3.05 

cm/sec 

Accelerometer Bias 

60 


Accelerometer Scale Factor 

500 

ppm 

Accelerometer Nonorthogonality 

10 

ppm 

Initial Gyro Angular Error 

42 

arcsec 

Gyro Bias Drift 

0.01 

deg/hr 

Gyro Scale Factor 

33 

ppm 

Gyro Nonorthogonality 

20 

ppm 

Gyro Random Walk Intensity 

0.025 

deg/hr 1 / 2 
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